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A method for illumination compensation in a system for processing of signals representing reflected light intensities from 
a surface includes the steps of receiving reflected light intensity information in the form of more than one intensity value (33) and 
receiving, associated Tvith each one of the reflected light intensity values, a value of at least one variable (37, 41). The method fur- 
ther includes remapping (42) each of the reflected light intensity values to one or more than one compensated reflected light in- 
tensity values, including remapping any of the reflected light intensity values equal to or less than a preselected threshold mini- 
mum reflected light intensity value associated with the value of each of the variables for that reflected light intensity value, to a 
minimum level output, and remapping any of the reflected light intensity values equal to or greater than a preselected threshold 
maximum reflected light intensity value, associated with the value of each of the variables for that reflected light intensity value, 
to a maximum level output 
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METHOD AND APPARATUS FOR COMPENSATION FOR 

NON-UNIFORM ILLUMTNATION 

Background of trhe Invention, 

This invention relates to compensation for non- 
5 uniformity in illumination of a surface in a system that 
takes information in optical form from the surface and 
converts such information to electronic form and in par«- 
ticular to compensating for variations in illumination in 
systems for reading labels on articles, such as packages 
10 or letters, being moved at high speed along a conveyor 
belt past a video camera. 

In the shipping of parcels and packages, it is 
necessary, at various points in the process, to obtain 
various information that may be encoded on a label. Such 
15 information includes address information, shipper informa- 
tion, and the like. Encoded labels also have application 
in other fields, such as manufacturing. Uses for informa- 
tion encoded on a IcUsel may include, by way of example, 
sorting packages further along the conveyor belt, deter- 
20 mining the quantity of items going to a particular loca- 
tion, and recording the movement of packages so as to 
provide updates in response to customer inquiries. The 
rapid obtaining of information from a label may be facili- 
tated by an automated system for converting optically 
25 encoded information on the labels to electrically encoded 
information. 

Xn modem terminals for sorting of packages, 
packages are loaded on a conveyor belt moving at speeds of 
up to 100 inches per second and moved to various points 
30 for sorting. Such conveyor belts are ordinarily three to 
four feet vide. A common width is 42 inches. Packages 
may be loaded onto such conveyor belts anywhere across 
their widths, and in any orientation. Packages of differ- 
ing sizes may be shipped along a single conveyor belt. As 
35 a result when a package passes through a zone for reading 
of label information, the label may be at any point across 
the width of the conveyor belt and at any height within a 
range of heights above the surface of the conveyor belt. 
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Reflected light intensity is comnonly converted 
to electronic form using a charge coupled device ("CCD"). 
An appropriate system of lenses focuses reflected light 
onto a CCD array. The CCD generates a signal, represent- 
ing an intensity value, on a pixel-by-pixel basis, for the 
zone for reading label information. Each pixel represents 
an area viewed by the CCD. Signal processing equipment 
then converts the intensity value for each pixel to a 
signal representing one of the possible optical properties 
of the label. Such optical properties may be vrtiite and 
black, as well as various shades of gray. The optical 
property is established by comparing the intensity value 
to one or more reference threshold values. 

It is difficult to provide an illumination 
15 system that can uniformly illuminate the entire zone 

through which labels may pass. As noted above, this zone 
is three to four feet wide and three to four feet high, 
m fact, the intensity of illumination may vary across the 
width of the zone at any given height. The intensity of 
20 illumination will ordinarily also vary with height. The 
pattern of intensity across the width of the zone is not 
the same at different heights. As a result, a CCD or 
other apparatus for converting reflected light intensity 
values into electrical signals may produce signals repre- 
25 senting different intensity values for the same optical 
property depending on the location of the label in either 
height or position across a conveyor belt. As a result of 
such variations in intensity values, actual threshold 
values for distinguishing different optical properties 
30 from one another will vary depending on the lopation of 
the label. If reference threshold values cannot be ad- 
justed to compensate for this effect, errors in reading 
labels will occur. For example, a white area of a label 
in a region of low illtamination may give a intensity 
35 signal that is below a fixed threshold for black or gray. 
This will result in errors in reading of labels. Such 
errors of course cause misreading of information on a 
label. 
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n%.-t^r^i>g »f t^he Invention, 

Accordingly, it is an object of this invention 
to provide a method and an apparatus for correcting sig- 
nals representative of information stored on a surface in- 
5 optical form detected under conditions of non-\inif ormity 
of illumination of the surface. 

It is a further object of this invention to 
provide such a method .and apparatus for use in a system in 
which labels are placed on packages of varying heights. 
10 It is a further object of this invention to 

provide a method and apparatus for correcting signals 
representing label information in which labels are placed 
on packages that are moving on a conveyor belt. 

It is a particular object of this invention to 
15 provide such a method and apparatus for use where the 
information is stored on the surface in white areas and 
black areas. 

Additional objects and advantages of the inven* 
tion will become evident from the detailed description of 
20 a preferred embodiment below. 
Summarv of the Invention, 

A method for illumination compensation for use 
in a system for reading surfaces having at least black and 
white areas, includes the steps of storing reference 
25 reflected light intensity information corresponding to at 
least black and white areas in the form of the plurality 
of intensity values, each one of the intensity values 
having an address defined by a value of at least one 
variable, calculating from the reference reflected light 
intensity information more than one black threshold value 
and more than one white threshold value each of which 
threshold values is associated with a selected value of 
each of the variables. The method further includes the 
steps of receiving opera *-tional reflected light intensity 
information in the form of more than one operational 
reflected light intensity value, each of the operational 
reflected light intensity values having an address defined 
by a value of each of the variables, remapping each of the 
operational reflected light intensity values to one of at 
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least two output compensated reflected light intensity 
values, including the steps of remapping any of the operar 
tional reflected light intensity values equal to or less 
than the black threshold value for its address to a black 
5 level output, and remapping any of the operational re- 
flected light intensity values equal to or greater than 
the white threshold value for its address to a white level 
output, and producing. an output signal including the 
output compensated reflected light intensity values. 

A method for illumination compensation in a 
system for processing of signals representing reflected 
light intensities from a surface includes the steps of 
receiving reflected light intensity information in the 
form of at least two intensity values, receiving, associ- 

15 ated with each of the reflected light intensity values, a 
value of at least one variable, and remapping each of the 
reflected light intensity values to one of at least two 
compensated reflected light intensity values, including 
the steps of remapping any of the reflected light intensi- 

20 ty values equal to or less than a preselected threshold 
minimum reflected light intensity value associated with 
the value of each of the variables for the reflected light 
intensity value to a minimxun level output and remapping 
any of the reflected light intensity values equal to or 

25 greater than a preselected threshold maximiim reflected 

light intensity value associated with the value of each of 
the variables for the reflected light intensity value to a 
maximum level output. 

An apparatus for illumination compensation in a 

30 system for processing of signals representing reflected 
light intensities from a surface includes means for re- 
ceiving reflected light intensity information in the form 
of more than one intensity value, means for receiving, 
associated with each one of the reflected light intensity 

35 values, a value of at least one variable, and memory means 
for remapping each of the reflected light intensity values 
to one of at least two compensated reflected light inten- 
sity values, including memory means for remapping any of 
the reflected light intensity values equal to or less than 
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a preselecrted threshold minimuia reflected light intensity 
value associated with the value of. each of the variables 
for the reflected light intensity value to a minimum level 
output and memory means for remapping any of the reflected 
5 light intensity values equal to or greater than a prese- 
lected threshold maximum reflected light intensity value 
associated with the value of each of the variables for the 
reflected light intensity value to a maximum level output. 
Brief Description of the Drawings, 

10 Figure 1 is a simplified perspective view of a 

system for use with a method and apparatus according to 
the invention. 

Figure 2 is a block diagram of a circuit accord- 
ing to the invention. 

15 Figure 3 is a partial view of an calibration bar 

for use in connection with the invention. 

Figure 4 is a flow chart of the overall process 
of calibration in connection with a method and apparatus 
according to the invention. 

20 Figure 5 is a graph showing reference x-eflected 

light intensity data versus horizontal position before 
simplification. 

Figure 6 is a flow chart showing the steps of 
reflected light intensity data simplification. 

25 Figure 7 is a graph showing reference reflected 

light intensity data versus horizontal position after data 
simplification. 

Figure 8 is a graph showing output reflected 
light intensity values versus input reflected light inten- 

30 sity values as remapped according to the invention. 
Detailed Description ef a Preferred Embodiment. 

With reference to Figure 1, there is shown an 
example of an apparatus for reading surfaces having infor- 
mation, such as reading label information on articles 

3*5 being moved on a conveyor belt. Such an apparatus may be 
used in connection with an illumination compensation 
appcuratus of the invention. There is shown a section of a 
conveyor belt 60. Conveyor belt 60 is moving in the 
direction of arrow 61. Box 62 is on the conveyor belt 60. 
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BOX 62 has height h. Box 62 has label 64 on a top surface 
thereof. For clarity, support members are not shown. It 
will be clear to those of skill in the art of design of 
mechanical supports that mechanical supports can easily be 
provided . 

In moving on conveyor 60 in the direction shown 
in the arrow, an article first encounters height sensing 
means 66. Height sensing means 66 includes two vertical 
apparatuses 68 and 70 disposed opposite each other across 
the belt. Apparatuses 68 and 70 may be any one of a 
number of commercially-available light sensors. One of 
apparatuses 68 is an emitter, and the other is a detector. 
Height sensing apparatus 66 is used to ascertain the 
height h of package 62. 
3^5 A package 62 next moves through illumination 

zone 72. Illumination zone 72 is a generally rectangular, 
planar region of high-intensity illumination provided 
across the width of the belt up to a selected height above 
the surface of the belt. The illumination is provided by 
20 illumination apparatus 74, which is located over the 

illumination zone. Light reflected from a box 62 proceeds 
upward to mirror 76 and then is reflected to light-detect- 
ing device 78. Vertical uprights are provided on each 
side of the illumination zone 72. Each vertical upright 
25 has a number of equally-spaced horizontally-protruding 

pegs 82 extending therefrom. Each peg 82 corresponds to a 
peg 82 at the same height on the opposite one of the 
vertical uprights 80. The pegs 82 support the calibration 
bar. The use of the calibration bar is explained below in 
30 connection with the discussion of calibration. 

The circuitry of the invention will now be 
explained with reference to Figure 2. Two types of infor- 
mation are received by the apparatus of the invention. A 
first type of information is reflected light intensity 
35 information from a light-detecting device. In a preferred 
embodiment, the light-detecting device is a CCD. This 
reflected light intensity information is transmitted along 
line 30 to initial gain unit 31. Reflected light intensi- 
ty data is transmitted by the CCD as signals corresponding 
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to discrete pixels. The reflected light intensity infor- 
mation is in the form of a series of discrete values or 
reflected light intensity values. Each discrete value 
corresponds to a particular position along a scan line. 

5 A second type of information is information 

relating to the values of the variables, or variable value 
information, upon which the illumination depends. In the 
disclosed embodiment of the invention, there are two 
variables upon which the illumination depends. One such 

10 variable is height. Height information is transmitted 
along line 34 to a computer 36. 

The second such variable is horizontal position. 
The horizontal position of any item of reflected light 
intensity is preferably derived from a clock signal asso- 

15 ciated with the reflected light intensity information 

transmitted by the light-detecting device. A clock signal 
from the light-detecting device is transmitted along line 
38 to position logic xinit 40. The light-detecting device 
also transmits a signal to the position logic unit 40 

20 representing the beginning of a new scan line. Prom the 
clock signal and the new scan line signal, position logic 
unit 40 provides a position signal. The unit may be an 
Altera 5032. The providing of such a signal is conven- 
tional in the art. The position signal is sent simulta- 

25 neously with each intensity value. The position signal 

counts from zero on each scan line across the illumination 
zone. The position signal value is changed a certain 
number of times in each scan line depending on detail of 
compensation in the apparatus. This is explained below in 

30 more detail where the look-up tables are discussed. 

Computer 36 is shown in Figure 2. As noted 
above, an input of computer 36 receives height informa- 
tion. Computer 36 has the function of processing height 
information in the system of the invention. One output of 

35 computer 36 is coupled to an input of initial gain unit 

31. This output is a gain signal. In a preferred embodi- 
ment, this gain signal is an 8-bit digital signal. This 
gain signal is used as an adjustment for gain variation 
relating only to height. It has been found that reflected 
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light intensity varies greatly witli height. It is be- 
lieved that this variation with height results from two 
causes. First, as the distance from the illumination 
means increases, the intensity of illumination falls off 
5 rapidly. Second, as the distance from the focusing means 
of the light detecting device increases, a given area on a 
surface corresponds to a smaller area on the light-detect- 
ing device. Accordingly, independent of intensity of 
illumination, the same optical property impinging on the 
10 CCD will cause fewer photons to strike the surface of the 
CCD for a package of lesser height than for a package of 

greater height. 

The output of computer 36 is a height signal. 
In a preferred embodiment, the height signal is a two-bit 

15 digital signal. Thus, the value of the height signal 

indicates four ranges of height in the illumination zone. 
The height signal output is coupled to the input of the 
illumination compensation circuitry 42. 

The output of initial gain adjustment 31 is a 

20 gain-adjusted analog reflected light intensity signal. 

The output of initial gain adjustment 31 is coupled to an 
input of analog-to-digital conversion circuitry 32. The 
output of analog-to-digital conversion circuitry 32 is a 
digital reflected light intensity signal. In a preferred 

25 embodiment, the digital reflected light intensity signal 
is an 8-bit signal. This output is coupled to an input of 
illumination compensation memory device 42. 

The output of the position logic 40 through line 
41 is coupled to an input of the illumination compensation 

30 memory device 42. This signal represents the variable 
horizontal position. In a preferred embodiment, this 
position logic output is a four-bit signal. 

It will thus be understood that the illumination 
compensation memory device receives reflected light inten- 

35- sity value information along line 33. This information is 
provided as a series of individual reflected light inten- 
sity values. Each reflected light intensity value has an 
address defined by variable value information for that 
reflected light intensity value. In Figure 2, the values 

-8- 
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of two va^riables define "the address of each reflected 
light intensity value. 

During calibration of an' apparatus according to 
the invention, the reflected light intensity value infor- 
5 mat ion represents light reflected from a surface having a 
known, or reference, optical property. Such reflected 
light intensity value information is referred to as refer- 
ence reflected light intensity value information. 

During operation of an apparatus according to 

10 the invention, the reflected light intensity value infor- 
mation represents light reflected. from a surface of un- 
known optical properties. Such reflected light intensity 
value infoanaatlon is referred to as operational reflected 
light intensity value information. 

15 The initial gain adjustment apparatus will now 

be explained in greater detail. An 8-bit digital height 
signal is provided from the computer. The signal is input 
to a multiplier such as a CA 3338. The output of this 
logic unit is a gain-adjusted signal containing reflected 

20 light intensity value information. This is an analog 

signal. The amplitude of the output signal is determined 
by the height signal provided by the computer. Xn the 
preferred embodiment, a chip is provided that has 256 
possible gain adjustments. 

25 The height signal is generated by the computer 

depending on the height of the package. In one embodi- 
ment, a height sensor is employed, over a range of 28 
inches, with a sensitivity of one-half inch. Thus, in 
that embodiment 56 different height readings are possible. 

30 In another embodiment, the height sensor is used over a 
range of 36 inches, thus making 72 different height read- 
ings possible. The computer has a look-up table stored in 
memory of gain signal values corresponding to each possi- 
ble height value. The entries in the look-up table are 
3*5 obtained during calibration of the apparatus. The cali- 
bration of the apparatus is explained in detail below. 

The output of the gain adjustment unit is cou- 
pled to an analog-to-digital conversion unit. The unit 
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may be, for example, a CA3 318. This cHip produces an 8- 
bit digital intensity signal • 

It will be seen that the illximination compensa- 
tion memory device has input lines dedicated to two types 
5 of input. A first type of input is uncompensated reflect- 
ed light intensity data. This data is transmitted along 
line 33 from A/D conversion unit 32. The uncompensated 
reflected light intensity data is in the form of a series 
of intensity values. 

Kie other type of input is variable value data. 
Each intensity value of uncompensated reflected light 
intensity data has an address defined by a value of each 
variable. Two input lines for variables are provided, 
one input line is for height information. The other 

15 variable input line is for horizontal position informa- 
tion. Thus each input intensity value is associated with 
a particular height value and horizontal position value. 

The operation of the illtimination compensation 
memory device will now be explained in more detail. The 

20 memory device maintains a set of look-up tables. The 
memory device uses the tables to remap input intensity 
values to compensated intensity values. 

Each table corresponds to a particular combina- 
tion of the variables. In a preferred embodiment, each 

25 table corresponds to a particular range of height and a 
particular range of horizontal location. The number of 
tables is dependent on the degree of detail desired in 
compensating. In a preferred embodiment, the illumination 
zone is divided into 32 segments of equal size across the 

30 width of the belt and 4 zones of equal size from the level 
of the belt to the top of the illumination zone. Conse- 
quently, 128 zones are provided. Thus, 128 tables are 
provided for adjusting the received intensity values. The 
table which is to be used for a given reading is deter- 

35 mined by the height value and position value for the 

reading. Each table has a six-bit address. The height 
value is set. for each entire label. The height signal 
provided to the illumination compensation memory device is 
accordingly a two-bit signal. The position value only 

-10- 
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changes 16 times in every sweep of the CCD across the 
width of the belt. The position signal is provided to at 
least two memory devices corresponding to left and right 
sides of the illxamination zone. In a preferred embodi- 
5 ment, there are four channels provided. Two channels 

receive reflected light intensity value information corre- 
sponding to the left half of the illumination zone. Two 
other channels receive reflected light intensity value 
information corresponding to the right side of the illumi- 

10 nation zone. Thus the position value has 16 possible 
values, while 32 horizontal segments are provided. The 
same position signal corresponds to two position variable 
values, depending on which channel received the signal. 
The position value signal is thus a four-bit signal for 

15 each half of the illumination zone. Each time the posi- 
tion value changes, a new table is designated on each half 
of the illumination zone. 

In the preferred embodiment, a CCD imager is em- 
ployed having 4,096 pixels in a scan line. Thus each 

20 table remaps 128 reflected light intensity values. 

The operation of the look-up tables in the 
memory devices is best understood with reference to Figure 
8. Figure 8 is a graph of output versus input for three 
horizontal segments at a single height range. Figure 8 

25 thus represents remapping of input to output for three of 
the 32 tables in a single height range. 

Line 802, which is a staircase curve, represents 
the remapping of the input values to the output values by 
one table. Line 802 corresponds to segment 2, which 

30 remaps values of pixels 128-255. All input digital inten- 
sity values less than a black level threshold value are 
remapped to a black output level. In the illustration of 
Figure 8, the black output level is zero. The black level 
threshold value for line 802 is approximately 102. All 

35 input values above a white threshold are remapped to a 
white output level. For the segments illustrated in 
Figure 8, the white output level is the maximum output 
level of the memory device. The white threshold for line 
802 is approximately 192. 
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Intensity values above the black level threshold 
and below the white level threshold are typically remapped 
according to a quantized linear function. Accordingly, an 
intensity value is remapped as follows. The output inten- 
5 sity value is equal to the sum of the black output level 
and the product of the absolute value of the difference 
between the black output level and the white output level 
by the quotient of the absolute value of the difference 
between the input value and the black level threshold 

10 divided by the absolute value of the difference between 
the white level threshold and the black level threshold. 
Remapping of the intensity values in this manner provides 
greater dynamic range for the data. Input values from a 
relatively compressed range between the black level thres- 

15 hold and the white level threshold are remapped to a wider 
range between the black level output and the white level 
output. For example, for segment 2, which corresponds to 
line 802, input intensity values from 105 through about 
190, are remapped to values from 15 to 224. 

20 The look-up tables for segments 8 and 14 are 

similar. The values of the black level threshold and the 
white level threshold are different for each segment. For 
example, for segment 9, represented by line 804, the black 
threshold value is about 118, and the white threshold 

25 value is about 220. For segment 14, represented by line 
806, the black threshold value is about 126, and the white 
threshold value is about 230. 

It should be noted that the remapping method of 
the invention is not confined to use with white and black 

30 threshold values. In fact, black is merely an example of 
a minimum reflected light intensity value. It is antici- 
pated that the minimxam reflected light intensity signal 
received by the memory device of the invention will corre- 
spond to a black surface. This is because it is antici- 

-95 pated that the system of the invention will be used with 
packages having labels printed with black areas. It is • 
possible, in a different use of the invention, that a 
minimum reflected light intensity value would correspond 
to an optical property of the surface other than black. 
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In such event, a reference reflected light intensity value 
would be obtained using a calibration bar having areas 
with whatever optical property represented the minimum re- 
flected light intensity value. This reference reflected 
5 light intensity value would be designated a threshold 

minimum reflected light intensity value. The memory means 
would store a look-up table designed to remap any opera- 
tional reflected light intensity value equal to or less 
than the threshold minimum reflected light intensity value 
10 obtained from the reference to a minimvim output level. 

Similarly^ it is possible to use the invention 
in a context in which there is a maximum reflective light 
intensity value that corresponds to an optical property of 
a surface other than white. In such a system, there would 
15 be a surface used for calibration that had the optical 
property corresponding to the maximiun reflected light 
intensity value used for determining a reference maximum 
reflective light intensity value. This reference reflect- 
ed light intensity value would be designated a threshold 
20 maximum reflected light intensity value. The memory 

device according to the invention would be programmed with 
a look-up table that remapped all operational reflected 
light intensity values equal to or greater than the thres- 
hold maximum reflected light intensity value to a maximum 
25 output level. 

It will be seen that in the general case, any 
operational reflected light intensity value falling be- 
tween the minimtim reflected light intensity value thres- 
hold and the maximtim light intensity value threshold would 
30 be remapped to an output value between a minimum output 

value and a maximtim output value. For example, the values 
of operational reflected light intensity value falling 
between the minimim reflected light intensity value thres- 
hold and the maximum reflected light intensity value 
-35 threshold can be remapped, to an output value between a 

minimum output value and a maximiim output value. However, 
the intermediate values could be remapped according to any 
desired function. 
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AS shown in Figure 2, the output values are 
transmitted by the meiaory device along an output line. 
The output signal is a series of compensated or remapped, 
reflected light intensity values. These nay be four-bit 
5 values, with the four higher bits of an eight-bit output 
line used. 

The process of calibration of the illumination 
compensation apparatus will now be described. It is 
necessary to carry out these calibration steps whenever 

10 changes are made that may affect the illumination. For 
example, replacement of a bulb or adjustment of a light 
requires that the illumination compensation apparatus be 
recalibrated. Recalibration includes reprogramming the 
memory device with tables for illumination compensation by 

15 remapping reflected light intensity values. 

An overview of the calibration process may be 
understood with reference to Figure 4. The first step is 
adjustment of the attenuation gain, as shown by box 502 in 
Figure 4. The next step is the step of obtaining refer- 

20 ence reflected light intensity values or profiles, shown 
by box 504 in Figure 4. The pattern of reflected light 
intensity values of any uniform optical property across 
the width of the zone is referred to as a "profile." The 
next step is the optional step of a hardware failure 

25 check, shown in box 506 in Figure 4. The next step is the 
step, recpiired in a preferred embodiment of the invention, 
of combining certain reference reflected light intensity 
values, or reducing the number of profiles, shown in box 
508 in Figure 4. The next step is the step of simplifying 

30 the reference reflected light intensity value profiles, 
shown in box 510 in Figure 4. The next step is the step 
of dividing the reference reflected light intensity values 
or profile data into segments, shown in box 512 on Figure 
4. The next step is the step of finding the inverse of 
-35 the reference reflected light intensity values or profile 
values, shown in box 514 in Figure 4. The next step is 
the step of multiplying the simplified and inverse simpli- 
fied values to calculate a flat response, shown in box 
515. The next step is the step of generating tables for 
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s'torage in a nemory device. This is shown in box 516 in 
Figiire 4. The last step is the step, shown in box 518 of 
Figure 4, of storing the tables in the memory devices. 

In making these calibration adjustments, it is 
5 desirable to use a calibration bar 100 as shown in Figure 
3. Alternating white strips 102 and black strips 104 are 
provided on the calibration bar 100 perpendicular to the 
length of the bar. For example, each strip may be 0.027 
inches in width. The calibration bar 100 is sufficiently 

10 long that it can rest on pegs 82 of Figure 1, and extends 
across the entire belt. It will be seen that for obtain- 
ing the reference reflected light intensity values, a 
surface or surfaces must be provided having both the 
optical property corresponding to the minimim reflected 

15 light intensity value and the optical property correspond- 
ing to the maximtim reflected light intensity value. 

In the first step, shown as box 502, the attenu- 
ation gain is adjusted. The adjustments are made using an 
oscilloscope to monitor the levels. The goal in adjusting 

20 the attenuation gain is to assure that the black level is 
not below zero, and that the white levels are not saturat- 
ed. The calibration bar is placed at various selected 
heights. The step of adjusting the gain is repeated at 
each selected height. The gain value for each selected 

25 height is stored in the computer for later processing. In 
a preferred embodiment, the selected heights are 1 inch, 4 
inches, 8 inches, 12 inches, 16 inches, 20 inches, 24 
inches and 28 inches above the conveyor belt. In a system 
for reading surfaces having heights up to 36 inches, 

30 additional readings would be taken. 

The gain values set in the calibration process 
represent eight height points. In fact, in a preferred 
embodiment of the invention, there is provided height 
sensing means having the capability of detecting half-inch 
^5 height intervals. The points intermediate the eight 

selected heights are obtained by standard curve-fitting 
techniques. In a preferred embodiment, the identical 
curve-fitting technique and software that is used as de- 
scribed below in the step of simplification of the com- 
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bined raw reflected light intensity data, or combined raw 
profiles, is employed. It may be seen, in the appended 
computer program listings, that the listing for this 
curve-fitting technique includes lines of code devoted to 
5 dividing by four the gain values obtained so as to enable 
use of the same code for this curve- fitting operation. 

The next step in calibration is the collection 
of data for forming the information required to accomplish 
illumination compensation. The collection of profile 
10 data is preferably accomplished simultaneously with the 
step of gain adjustment. The readings given when the 
calibration bar 100 is placed across the illumination zone 

are recorded. 

An image of the calibration bar is created at 
each selected height. This image is raw, or uncompensat- 
ed, data. The image of the calibration bar provides a 
reference image for the black level and the white level at 
a large number of horizontal positions. This raw image or 
raw profile, or reference reflected light intensity value 
information, is transmitted to the computer. In the step 
of transmitting data to the computer, the hardware is 
configured so that the raw data is transmitted to the 
computer without compensation. In a preferred embodiment, 
each illumination compensation memory device is programmed 
25 during calibration so that the data output by the illumi- 
nation compensation memory devices is the same as the 
input. This is accomplished by removal of the memory 
devices, or PROM's, used in the prior run of the system 
with so-called "flat PROM's." These, "flat PROM's" do not 
30 adjust the data that is input. The output can be truncat- 
ed before being stored. It may be truncated by a Pixel 
Bus (Trademark of Analogic, Peabody, Massachusetts) inter- 
face board, which is part of an embodiment of a computer 
system which receives data from the PROM's through a 
-35 buffer, and transmits data in truncated form to the re- 
mainder of the computer system. 

AS may best be seen with reference to Figure 5, 
the profile when viewing a calibration bar, or reference 
reflected light intensity value information, is substan- 
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tially in the form of two lines. With reference to Figure 
5, the horizontal axis represents horizontal position 
across the illumination zone. The vertical axis represents 
reflected light intensity. In Figure 5, thirty-two data 
5 points, or pixels , have been chosen for display in the 
figure* The intensity values shown have been measured by 
placing an input of an oscilloscope on an intensity output 
line from the light-detecting device. The intensity value 
numbers are arbitrary. 

10 The line designated by 602 in Figure 5 shows the 

white levels. As may be seen, there is substantial varia- 
tion in the white levels. In particular, the white levels 
drop on the edges of the illvimination zone. The line 
designated by 604 represents the black levels. As may be 

15 seen, there is also considerable variation in the black 
levels. A profile similar to Figure 5 is obtained for 
each height interval at which the calibration bar is 
placed. The set of raw profiles are stored for further 
processing. Each raw profile is stored in a separate file 

20 in the computer. 

After collection and storing of the reference 
intensity values, there is an. optional step of checking 
the reference intensity value information data for evi- 
dence of hardware failure. A hardware failure in the 

25 system would make reading of labels with a reasonable 
degree of accuracy impossible. Such hardware failures 
include a failure of a portion of the illimination appara- 
tus, such as a bulb, or of a portion of the light sensing 
means, such as a failure of a CCD array. 

30 Checking the profile data for evidence of hard- 

ware failures is accomplished by computing a variance and 
standard deviation. The equation for the variance, var 
(Xx..n is the following: ^ a/ z 

Where X is the intensity value and N is the number of 
35 pixels in a profile. The standard deviation <r(^|..n) 
is defined by the equation of ^{%^..^)- /vAr(^c,..n) 
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m a preferred embodiment, a routine to calculate the 
variance and standard deviation is taken from W. Press, et 
al., "Numerical Recipes in C", (Cambridge university Press 
1988) at page 475. The listing for the software routine 
5 called "moment. c" is shown in the appended source code 
listings. 

The next step after the optional step of check- 
ing the raw profile data for hardware errors is the reduc- 
tion of the number of profiles obtained in the data col- 

10 lection step to the number of profiles actually used 

during illumination compensation. For example, in a pre- 
ferred embodiment, in the data collection step, eight 
profiles are collected and stored. In a preferred embodi- 
ment, the illumination zone is divided into only four 

IS vertical divisions. Consequently, the eight profiles must 
be combined into four. It has been found desirable to 
take readings at levels 1 inch, 4 inches, 8 inches, 12 
inches, 16 inches, 20 inches, 24 inches and 28 inches 
above the belt. It has been found desirable to use the 

20 highest and lowest profiles individually. This is because 
these profiles are the extremes of lighting, and usually 
quite different from an average of intermediate height 
profiles. In this embodiment, the 1 inch and 28 inch 
profiles are used without adjustment. The 4 inch, 8 inch 

25 and 12 inch profiles are averaged on a pixel-by-pixel 

basis. The 16 inch, 20 inch and 24 inch profiles are also 
averaged on a pixel -by-pixel basis. The result is four 
combined profiles. This is referred to as the combined 
reference reflected light intensity value information. 

30 The next step is the step of simplifying the 

reference intensity value data. This is accomplished in a 
preferred embodiment by applying curve-fitting techniques 
to the combined raw profiles. A variety of curve-fitting 
algorithms may be used. In a preferred embodiment, a 
-J5 second-degree polynomial curve-fit algorithm is used. The 
procedure to convert raw profile data is as follows. The 
data is represented on a graph. The y-axis of the graph 
corresponds to reflected light intensity value and the x- 
axis corresponds to a horizontal position in the illumina- 
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tion zone. The equation for finding "the fitted curve from 
the raw profiles is the following: 



m 
























1-1 
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The source code for this is shown in CURVE-FIT 
C, which is included herewith in the appended software 
5 listinigs. 

Referring to Figure 6, there is shown a flow 
chart depicting the steps of the process of data simplifi- 
cation employing a preferred curve-fitting technique. In 
the step illustrated by the box labeled 702 on Figure 6, 

10 there is shown the step of allocating arrays. The alloca- 
tion of arrays is a conventional preliminary step to 
performing vector and matrix mathematics. The next step 
is the step of inserting the X and Y values in the matri- 
ces and arrays. The X values are the pixel numbers, which 

15 correspond to horizontal position in the illiimination 

zone. The Y values are reflected light intensity values. 
This step is shown at box 704 in Figure 6. The next step, 
shown in box 706 in Figure 6, is the step of computing the 
X summations. The x summations are elements of the left- 

20 hand matrix in the above equation for the curve-fit. 

This step requires computation of sizmmations of an array 
of X-values raised to various powers. A preferred soft- 
ware subroutine for this computation is shown under the 
file name "sum_a*' in the appended source code listings. 

25 The next step in the process is the step, shown in box 708 
in Figure 6, of computing the XY summations. This step 
requires taking two arrays of data, the array of X-values 
and the array of Y-values, and calculating a summation. 
The values in the array of Y-values are raised to increas- 

30 ing powers. A software siibroutine for this step is shown 
under the file name "sum_ab" in the appended source code 
listings. The next step, shown at box 710, is transfer- 
ring the values obtained from the step of computing the X 
summations to a matrix. This is the matrix shown on the 

35 left hand side of the above ecpiation for* the second-degree 
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polynomial curve-fit. As may be seen from inspection of 
the above equation, it is necessary to rearrange the 
equation to solve for the vector (b,, b,, b,) . Thus, it is 
required to find the inverse of the matrix of X summa- 
5 tions. This step is performed in a preferred embodiment 
by first finding the inverse of the matrix of X summations 
on a column-by-column basis. This is done by computing 
the I.U decomposition of a row-wise permutation of the 
matrix. A preferred function for finding the LU decompo- 

10 sition is found in the above-mentioned "Numerical Recipes 
in C" at Chapter 2, page 43. The source code for a soft- 
ware routine is listed under file name "ludcmp.c" in the 
appended source code listings. This step is shown at box 
712 in Figure 6. With reference to Figure 6, it will be 

15 seen that the next step is solving the resulting linear 
equations shown in box 714. A function for solving such 
linear equations is disclosed in the above-mentioned 
"Numerical Recipes in C" , chapter 2, page 22. The source 
code for a software routine for this function is listed 

20 under the filename "lubJcsb.c" in the appended sotirce code 
listings. The next step after solution of the linear 
equations is calculating the coefficients of B, as shown 
in box 716 in Figure 6. once this step is complete, a set 
of appropriate cujrves has been derived. 

25 The result of the step of simplifying data 

using a second-degree polynomial curve-fit is illustrated 
in Figure 7. The graph shown in Figure 7 is similar to 
the graph shown in Figure 5. Thus, the horizontal axis of 
Figure 7 represents the position horizontally across the 

30 illumination zone. The vertical axis represents a value 
of intensity of reflected light. Line 802 near the top of 
the graph represents the reflected light intensity data, 
after the step of data simplification, for the white 
strips on the calibration bar. Line 804, extending across 

35 near the bottom of the graph, represents reflected light 
intensity values, after data simplification, for the black 
strips on the calibration bar. 
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After the step of processing the raw profile 
data into curve-fitted profile data, there are four sets 
of profile data. Each set of profile data represents a 
profile corresponding to one of four ranges of height 
5 values in the illumination zone. Each set of profile data 
is next broken into 32 segments as shown by box 512 in 
Figure 4* The 32 segments represent the 32 horizontal 
segments of the illumination zone. Each of the 32 seg- 
ments corresponds to 128 pixels across the illumination 
10 zone. 

The simplified reference reflected light inten- 
sity values, or simplified profile data, as shown by box 
514 in Figure 4, is then mapped to an inverse curve. In 
obtaining the inverse, each line is first rescaled by 
15 conventional techniques. 

The values contained in the inverse curve are 
used in programming the illumination compensation memory 
devices. The values must be placed in a format that oan 
be accepted by the memory devices. In a preferred embodi- 
20 ment, an Intel HCS**86 Hex file is created for each of the 
PROM'S. The data stored is in the form of a table having 
256 elements. 

The data is then converted to a form which is 
compatible with the various channels into which the data 
25 is divided by the CCD units. For example, in a preferred 
embodiment, there are four channels. These four channels 
represent respectively odd and even pixels on the left and 
right sides of the belt. Accordingly, the data generated 
is divided into appropriate segments for each channel. 
30 Once the data is divided in this manner, the data is 
stored in computer memory in four parts. The data is 
transferred from the computer to the memory device. In a 
preferred embodiment, the memory device is an ultra-violet 
light erassible programmable read-only memory (PROM) . A 
35 Cypress 7C254 may be used, for example. In order to 

transfer data from the computer to such a device, a PROM 
burner is provided. The input of the PROM burner is 
coupled to a port of the computer so that data may be 
transferred from the computer to the PROM burner. In 



wo 93/04442 



PCr/US92/07a54 



accordance with usual techniques, the PROM burner is 
connected to the PROM so as to bum the data into the 
PROM, once the data has been stored in the memory device, 
the memory device is then inserted into the appropriate 
5 board. The process of burning and inserting the PROM is 
repeated for each channel. When the memory devices have 
been programmed and are in place in the circuit, the 
apparatus of the invention may be operated. 

The look-up tfiJDles may also be maintained in a 

10 random access memory. Such a memory would be programmable 
from a system computer. 

It will be appreciated that there are consider- 
able variations that can be accomplished in a method and 
apparatus according to the invention without departing 

15 from its scope. As a result, although a preferred embodi- 
ment of a method and apparatus of the invention has been 
described above, it is emphasized that the invention is 
not limited to a preferred embodiment and that there exist 
other alternative embodiments that are fully encompassed 

20 within the invention's scope, which is intended to be 
limited only by the scope of the appended claims. 
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* File Nane: noaent.c 
* 

* File Function: Routine ^o compute ^e average, standard deviation 

* and varaiance of "the dat;a passed. 

* Notes: Function was taken fron Kumerical Recipe in *C' see * 

chapter 13 pg 475 

* this function is a subset of the s&onent function. 
* 

* History: 1986 Created 
* 

/include <iDath.h> 



void aoment ( data , m , n , index , ave , sdev , svar ) 
int a, n, index; 

float **data , »ave , *8dev , *svar ; 
{ 

int j; 

float s; 

void nrerror ( ) ; 



if (n 1) 
nrerror { "n must be at least 2 in MOMENT" ) ; 

ssQ . 0 ; 

for (j = l; j <= n;j-M-) 
s ^ data [ index] [j] ; 

*avees/n; 

♦svar = 0.0; 

for(j c= 1; j <« n; 

*svar fabs(s = da ta[ index] [j] - (*ave)); 



*svar /= n; 
*sdevssqrt(*6var) ; 
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* 

* fCopyright 1991 United Parcel Service of America, Inc. 

* All Rights Reserved. 
* 

* Confidential, Unpublished Property of United Parcel Service 

* Use and Distribution Limited Solely to Authorized Personnel. 
* 

* The use, disclosure, reproduction, modification, transfer, 

* or transmittal of this work for any purpose in any form or 

* by any means without the written permission of United Parcel 

* Service is strictly prohibited. 
* 



* File Name: cxirvefit.c 
* 

* File Function: To take the profile or sa gain values and apply 

* a second degree curve fit to the data. The 

* function is controlled by the variable type 

* which selects arrays size for the data inputs 

* again profile data or sa gain data. 

* The output of the function is the three 

* coefficients that best fit the data. 
* 

* Notes: 
* 

* History: JSK 4-June-91 Created 

* JSK 20-June-91 Fix x index for SA data 
* 

/include <stdio.h> 
/include <math.h> 
/include "nr.h" 
/include "nrutil.h" 

/define VECMAX 32 
/define ARRAKAX 3 
/define SUKAX 4 

float sum_a ( ) ; 
float sum^ab ( } ; 
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curvef it (data , type , index , b) 
float **data,*b; 
int index; 
char type; 

{ 

float *x , *y , **a , *xsum , ♦xy sum , *col , **ai , d , power ; 
int io f )Cf»r *indx; 
FILE ^fp; 

X = vector (l.VECMAX) ; /* allocate array space */ 

y = vector (l.VECMAX) ; 

xsum = vector C 1 , VECMAX) ; 

xysum = vector (1, VECMAX ) ; 

a = matrix(l,ARRAMAX,l,ARRAMAX) ; 

ai B 3iiatrix(l,ARRAMAX,l,AKKAMAX) ; 

col ^ vector (ir VECMAX) ; 

indx = ivector (1, VECMAX ) ; 



if (type »P") 
in = VECMAX; 

else 
n » VECMAX/4; 



/* Dealing with profile gains */ 
/* Dealing with SA gains */ 



if (type — •?•) /* Profile Gains */ 

for(i « 1; i <= n; i-M-) /* Get x values */ 
x[i) = i; 

} 

else /* Sa gains */ 

{ 

for(i « 1, j l; i <s= m; i-j-+, j 4-= 4) 
^xti]=j; /* Note Sa gains are increments 4 */ 

for(i = l; i <= a; i++) /* Get y values */ 

y[i] = data[index][i]; 

power =1.0; /* Compute the x summations */ 

for(i = l; i <= SXJMAX; i+-K) 
xsiam[i] = suin_a (x,m,power-*H-) ; 



power =1,0; /* compute the xy summations */ 

xysum[l] = sim_a (y,m, power) ; 

for(i « 2; i < SUMAX; i++) 
xysum[i] = sum_ab(x,y,m,power++) ; 
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a[l][l] 
a[X][2] 
a[l]C33 
a[23[l] 
a[23[2] 
a[2][3] 
a[3][l] 
a[3][2] 
a[3][3] 



(float) 

X5U2a[l] 

xsim[2] 
xstm[l] 
xsiaD[2] 
xsim[3] 

XSX23n[2] 

X5i2in[3] 

XS122tt[4] 



a; 



/* Fill in lAe A matrix */ 



ludcinp(a,AMlAMAX, indXr &d) ; 
for(j = l; j <= ARRAMAX; jn 

^ for{i = l; i ARRAMAX; 

colti] = 0*0; 
coltj] = 1.0; 

lubksb ( a , ARRAMAX , indx , col ) 
for{i = l; i <« ARRAMAX; 
ai[i][j] = col[i]; 

} 



/* find inverse column by column 
/* perform LU decomposition */ 
)/* Math functions ludcmpO & */ 
lubksbO obtain from */ 
/* Numerical Recipes */ 



i++) 



/* 
/* 



Solve 
A x X 



n linear 
= B */ 



ec[uations 



for(i = 1; i <= ARRAMAX; i-M-) /* Now find coefficients of b */ 
{ 

bti] = O.O; 

for(j l; j <= ARRAMAX; j++) /* b = A xnverse * y */ 
b[i] ai[i][j] * xysum[j]; 



free_i vector ( indx , 1 , VECMAX ) ; 
free vector(col,l, VECMAX) ; 
free"matrix(ai,l, ARRAMAX,!. ARRAMAX) ; 
f ree^matrix( a, 1, ARRAMAX,!, ARRAMAX) ; 
f ree_vector (xysum^ l , VECMAX) ; 
f ree_vector (xsum , 1 , VECMAX) ; 
f r ee_vect or (y,l, VECMAX) ; 
free vector (x, 1, VECMAX ) ; 



/* Deallocate space */ 
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* File Name: stun a 

* File Function: To take an array of data and find a susunation * 

raised to a variable power 

* Notes: 

* History: JSK 4-June-91 Created 

float s\2in_a(a,lifflit,power) 
float *a, power; 
int limit; 

{ 

int i; 

float s\xm = 0.0; 

for(i =1; i <= limit; i++) 
sum +« pow(a[i] , power) ; 

return sum; 

} 

/***************************************************^************ 

* File Name: sum ab 

* File Function: To take two arrays of data and find a sximmation 

* raised to variable power 

* Notes: 
* 

* Historj': JSK 4-June-9i Created 
* 

float sum_ab (a /b, limit, power) 
float *a7*b, power; 
int limit; 

i 

int i; 

float sum = 0.0; 

for(i = 1; i <= limit; i++) 
sum (pow(a [i] #power) * b[i]); 

return sum; 

} 
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* File Name: lubksb.c 
* 

* File Function: Solves the set of n linear equations A*X = b 

* Notes: Function was taken from Numerical Recipe in 'C* see • 

* chapter 2 pg 44 

* History: 1988 Created 
* 



void lub)csb(a,n,indXrb) 
float ♦*a,bt]; 
int n,*indx; 
{ 

int i,ii=o,ip,D; 
float sum; 

for(i = l;i <= n;i++) 
{ 

ip = indx [ i ] ; 
sum - b[ip] ; 
b[ip] « b[i]; 
if (ii) 

for(j = ii;j <= i-l;D++) 
SUB -= a[i] [j]*btj] ; 
else 
if (Sim) 
ii = i; 
b[i] = sum; 

} 

for(i = n;i >= l;i — ) 
{ 

sum = b [ i ] ; 

for{j - i-i'l;j <= n;j-h+) 

sum -= a[i][j]*b[j]; 
b[i] = Bum/a[i] [i] ; 

} 
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* File Naine: ludcmp.c 

* File Function: Given an nxn matrix routine replaces it by the LU 

* decomposition of a row vise permutation of * 
itself. 

* 

* Notes: Function was taken from Numerical Recipe in 'C" see * 

chapter 2 pg 4 3 

* History: 1988 Created 

/include <math.h> 

/define TINY l.Oe-20; 

void ludcmp(a,n,indx,d) 
inr n,*indx; 
float **a,*d; 
{ 

int i,imax,j,k; 

float big,dumr sum, temp; 

float *w,*vector () ; 

void nrerrorO ,free_yector () ; 

w^vector (i,n) ; 
*d«1.0; 

for {i«l;i<«sn;i++) { 
bigso.o; 
for j<»n; 

if ( (temps^f abs(ati] [ j]) ) > big) big=temp; 
if (big == 0.0) nrerror ("Singular matrix in routine 
LTOCMP»») ; wti]=1.0/big; 
} 

for (3=1; j<=n;j-i-+) { 

for (i«l;i<j;i++) { 
sum»a[i] [j] ; 

for (k«=l;k<i;k+-f) siim -« a[i] [k]*atk] [ j] ; 
a[i] [j3=sum; 

} 

big>=0»0; 

for (i=j;i<=n;i-M.) { 
sum=a[i3 [ j] ; 
for (k=l;k<j ;k+-f ) 

sum a[i][k]*a[k][j]; 
ati] [j)=Bum; 

if ( (dum=w[i]*fabs(sum) ) >= big) { 
bigscium ; , 
imax=i ; 

} 

} 
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if (j !- imax) { 

for {lc=l;)c<=n;k++) { 
dun=aCimax] [)c) ; 
a[iiBax] [)c]=a[ j] [Jc] ; 

} 

*d = -(*ci) ; 
wtimax)=w[ j] ; 

} 

indx[ jj^imax; 

if (atjJtj] = 0-0) a[3]C3]=TINY; 
if (j 1= n) { 

duin-1.0/(atD3 [D]) ; 

for (i=j+l;i<=n;i++) a[i][j] *- dum; 

} 

} 

f r ee_vector ( w , 1 , n ) ; 

} 



#undef TINY 
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/* 

* File : genprom . c 

* Aut:h: Robert K Pekarek 

* Date: 03 Jan 1990 
* 

* Copyright (c) 1990 United Parcel Service. 

* Purpose: Generate prom object files for SA correction tables. 
* 

* Edit: Jim S Kunicki (JSK) added /define min marco to compile 

* under Unix 

* Changed progan to be a callable function 
*/ 

#def ine HALVES 2 
/define FIRST 0 
/define SECOND 1 
/define HZ 4 
/define PZ 16 
/define DATABITS 256 
/define SIZE 16 
/define VERSION 1 
/define REVISION 0 

/define min (a, b) (((a) < (b) ) ? (a) : (b) ) 



/include <6tdlib.h> 
/include <stdio.h> 
/include <string.h> 
/include <ctype.h> 
/include <math.h> 
/ifdef KSDOS 
/include <process.h> 
/endif 

/*«**********************»***.**^***********************^***^*^^ 

extern void main(int argc, char **argv) ; 

static FILE *GetInputFile(int argCount, char **stringVector) ; 
static void GetFSOpsRND(riLE *fp, int *full, float ♦preScale, int 
♦quantize, 

int ♦inputScale) ; 
static void BuildRawOutOff set (FILE *fp, float *offl, float ♦rawl, 

float *of f 2 , 

float *rav2) ; 

static void CalcProm( float *offl, float *rawl, 

int full, float preScale, int quantize, 

int inputScale, 

unsigned char ^result, float *scale) ; 
static void WritelntelHex (unsigned char ^result, int aLen, char 
*outFileName) ; 

static void WriteText (unsigned char *result, char *outFileName) ; 
static void WriteScaleFactors (float *scale, int aLen, char 
*outFileName) ; 
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static FILE *Ge1:Input:File ( ) ; 
static void BuildRawOutOf f set ( ) ; 
static void CalcProm(); 
static void WritelntelHex ( ) ; 
static void WriteText ( ) ; 
static void WriteScaleFactors( ) ; 
static void GetFSOpsRND( ).; 

/* Global arrays */ 

float rawOut[HALVES][HZ][PZ], of f set [HALVES] [HZ] [PZ] ; 

float prescale, scale[HALVES] [HZ] [PZ] ; 

unsigned char r esults [ HALVES] [HZ] [PZ].[DATABITS] ; 



* main - get info and build prom files 
♦/ 

genpron ( ) 

{ 

FILE *fp; 

int f ullScaleOutput , roundToNear, inpScale; 

fprintf (stderr, "Genprom version *li.%li\n", VERSION, 
REVISION) ; 

/* open file for input ♦/ 
if ( (fp = GetlnputFileO ) = NDLL) 

^ fprintf (stderr, •'Usage: <reguired file naiae>\n"); 
return; 

} 

/* Get full scale output, pre-scale, and quantize factor */ 
GetFSOpsRND ( f p ^ & f ullScaleOutput , fipreScale , fcroundToNear , 
&inpScale) ; 

/* Build input offset and raw data arrays */ 

BuildRawOutOffset(fp,6offset[FIRST] [0] [0] , firawOut [FIRST] C01[0] , 

&of f set [ SECOND] [ 0 ] [ 0 ] , fcrawOut [ SECOND ] [ 0 ] [ 0 ] ) ; 

/* Close input file */ 
f close (fp) ; 

/* Calculate scale factors and results */ 
CalcProiD(&offset[ FIRST] [0] [0] , &rawOut[FIRST] [0] [0] , 

f UllScaleOutput, preScale, roundToNear, inpScale, 
fcresultsf FIRST] [O] [0] [0] , &scale[ FIRST] [0] [0]) ; 
CalcPrott(&offset[ SECOND] [0] [0] , firawOut[ SECOND] [0] [0] , 

f UllScaleOutput, preScale, roundToNear, inpScale, 
&results[ SECOND] [0] [0] [0] , &scale[ SECOND] [0] [0]) ; . 
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/* Write proa data to files */ 

WriteIntelHex(&resultst FIRST] [0] [0] [0] , P2*HZ*DATABITS, 
"AB.PRM") ; 

WriteIntelHex(&results[ SECOND 3 [0] [0] [0] , P2*HZ*DATABITS, 
"CD.PRM") ; 

/* Write prom data to files */ 

WriteText(&results[ FIRST] [0] [0] [0] , ••AB.TXT'*) ; 
WriteText(&results[ SECOND] [0] [0] [0] , "CD. TXT") ; 

/* Write out scale factors */ 

WriteScaleFactors(&scale[ FIRST] [0] [0] , P2*H2*HALVES , 
"SF.TXT"); 

/* Acknowledge completion */ 

fprintf (stderr, "Done generating prom files and scale 
f actors \n"); 
} 

/* 

* GetlnputFile - verify argiment count and retrieve file name 
*/ 

/* */ 

static FILE *GetlnputFile ( ) 
{ 

FILE *tf ; 

/* Return invalid argument count */ 



if (argcount 1= 2) 

{ 

return NULL; 

} 

/* Open file note hard coded for now */ 
if ((tf = fopen( "DATA. OUT", "rt")) -= NULL) 
{ 

perror ("") ; 

> 

return tf; 

) 

/* 

♦ GetFSOpsRND - read scalers from input file in designed format 
•/ 

/* */ 

#ifdef SPARC 

^define _MAX_PATH 144 /* Note the was a MICROSOFT "C" define 

within */ 

#endif /♦ stdlib.h, SUN "C" is not supported */ 
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static void GetFSOpsRND ( f p , full , preScale , quantize , inputScale) 
FILE *fp; 

int *f ull , *quanti2e , *inputScale ; 

floar •preScale; 

{ 

char buffer [_MAX_PATH ] ; 
int rc; 

/* Get preScale from first line */ 

if (f gets (buffer, _MAX_PATH, fp) = NULL) 

{ 

perror ("") ; 
eacit(l) ; 

if ((rc = sscanf (buffer, "PS = %f", preScale)) != 1) 

{ 

perror ("") ; 
exit(l) ; 

} 

/* Get full scale output from second line */ 
if (fgets( buffer, _MAX_PATH, fp) -= NULL) 
{ 

perror("") ; 
exit(l) ; 

if ((rc - BBcanf (buffer, "FSO = *i", full)) l= l) 
{ 

perror ("") ; 
exit(l) ; 

} 

/* Get round to nearest from third line */ 
if (fgets( buffer, _MAX_PATH, fp) = NULL) 
< 

perror ("") ; 
exit(l) ; 

if ((rc = sscanf (buffer, "RND = %i", quantize)) != 1) 
{ 

perror ("") ; 
exit(l) ; 

> 

/* Get input scaling from fourth line */ 
if (fgets( buffer, _MAX_PATH, fp) «=- NTJLL) 

perror ("•') ; 
exitd) ; 

if ((rc = sscanf (buffer, "IS = *i", inputScale)) != 1) 
{ 
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perror ; 
exit(l); 

} 

} 

/* 

♦ BuildRa%rt)utOff set: - read raw outputs and offsets from input file 
*/ 

/♦ */ 

static void BuildRawoutOffset (fp,off l,rawi,off 2,raw2) 
FILE *fp; 

float *of f 1 , *rawl , *of f 2 , *raw2 ; 

{ 

#def ine NUMPAEKS 2 

/define LONGBUFFER 300 

#define STARTPOS 23 

static char MatchString [ ] = "Diag"; 

char buffer [LONGBUFFER] , *rc, ♦bufferPtr; 

float *ptrl, *ptr2; 

int hz, pz, count, buf f erOf f set; 

/* Find first line containing the identifying string */ 
do 

{ 

rc = f gets (buffer, LONGBXJFFER, fp) ; 

} while ({rc != NULL) && strncmp (buffer, MatchString, 
strlen (MatchString) ) ) ; 

if (rc = KULL) 
{ 

perror("No data found"); 
exit(l); 

} 

/* For all height zones in the file, and not at end of file 

*/ 

for (hz = 0; (hz < HZ) (rc 1= ITOLL) ; ++hz) 
{ 

/* For both sets of input data */ 

for (ptrl offl, ptr2 « off 2 + P2 - i, count » O; 
count < KUMPARMS; 

ptrl « ravi, ptr2 «= raw2 + PZ - i, -M-count) 

{ 

/* Locate beginning of data on line */ 
bufferPtr = ^buffer [STARTPOS] ; 

/* Get the data for the left side */ 
for (pz = 0; pz < PZ; ++pz, ++ptrl) 
{ 

sscanf (buf f erPtr, " *f,%n", ptrl. 
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&^uf f erof fset ) 



&bufferOffset} 



} 



bufferPtr +- buff erOff set; 



/* Get the data for the right side ♦/ 
for (pz = 0; pz < PZ; ++pz, — ptr2) 

^ sscanf (buf f erPtr , " %f,%n", ptr2 , 

bufferPtr buf f erOf f set; 

} 

/* Get neart line */ 

rc = f gets (buffer, LONGBUTFER, f p) ; 



/* Increment data pointers */ 
offl ^-^ PZ; 
off 2 PZ; 
ravi PZ; 
rav2 PZ; 



/* 

* CalcProm - compute prom index values from given empirical data 
*/ _ 

i* t I *^tic void 
CalcProm ( of f l , rawl , full , preScale , quantize , inputScale , result , sea le ) 

float *of fl, *ravl, preScale, * scale; 

int full , quantize , inputScale ; 

unsigned char *result; 

{ 

int outerLoop, sel; 

float tempScale, tempStore; 

/* For each of the height and position zones */ 

for (outerl-oop =0; outerl*oop < H2 * PZ; -M-outerLoop) 

{ 

/* Compute and store the scale value */ 
tempScale » *scale-M- « 

(full / inputScale) / (preScale * *ravl - *of f 1) ; 

/* For each data selector ♦/ 

for (sel « 0; sel < DATABITS; ++sel) 

{ 

/♦ Compute the desired output and store to the 

result area */ . , 

tempStore = sel - *offl * inputScale; 

if (tempStore < (float) 0) tempStore = (float) 0; 

tempStore *= tempScale; 
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if (tempStore > (float) full) tempStore 

(float) full; 

tenpStore = (float) floor (tempStore / 

(float) quantize (float) 0.5) * 

(float) quantize; 

*result++ = (unsigned char) tempst ore; 

} 

/♦ Increment pointers */ 

ravl-H- ; 

offl++; 

} 

} . . 

/* 

* WritelntelHex - put character data out to given file in the 
intel hex 

* format 

~ * 

static void Kr itelntelHex (result , aLen , outFileName ) 
iinsigned char *result; 
int aLen; 

char *outFileName; 

{ 

FILE *fp; 

int count, tempLen, chkSum, dataCount; 
/* Open file */ 

fp « fopen( outFileName, **vt*"} ; 

/* For each data group of SIZE */ 
count » 0; 
while (aLen > 0) 

{ 

/* Put out intel Hex line ♦/ 
tempLen » min(SXZ£, aLen); 

chkSum « 0x7000 - tempLen • (count / DATABITS) - (count 
% DATABITS); 

fprintf(fp, ••:%O2X%O4X00»', tempLen, count); 

/* Put out data */ 

for (dataCount = 0; datacoiint < tempLen; -♦-+dataCount) 
{ 

fprintf(fp, "%02X", *result) ; 
ch)csum -« *result++; 

} 

/* Put out checksum */ 

fprintf(fp, ••%02X\n", chXSum % DATABITS); 

aLen tempLen; 
coiint tempLen; 
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} 



} 

/* Put out end record */ 
fprintf(fp, ":0000000lFF\n") ; 

fclose(fp) ; 



* WriteText - put character data out to given file in text format 
*/ 

/* */ 

static void WriteText (result rOUtFileName) 
unsigned char ^result; 
char *outFileNase; 

{ 

FII-E *fp; 

int hz, pz, datacoiint; 
/* Open file */ 

fp = fopen(outFileNaine, "vt") ; 

/* For each data group */ 
for (hz « 0; hz < HZ; -M-hz) { 
for (pz e 0; pz < PZ; HH- pz) { 

fprintf(fp, "\n HZ = %li PZ - %2i\n", hz, pz) ; 

/♦ Put out data */ ^ ^ « 

for (datacount « 0; dataCount < DATABITS; ++dataCount) 

^ fprintfCfp, -% 3u *result-w-) ; 

if (i ((datacount + 1) % SIZE)) fprxntf{fp, "\n") ; 

} 

}} 

f close (fp) ; 

} 

WriteScaleFactors - put scale factors to given file in ASCII 
format 
*/ 

/* */ 

static void WriteScaleFactors (scale, aLenrOUtFileName) 
float *scale; 
int aLen; 

char «outFileName; 

{ 

FILE ♦fp; 
int count; 

. /* Open file */ 
fp = fopen(outFileName, "wt") ; 
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fprintf(fp, "Channels A/B\n"); 

/* Write all scale factors */ 

for (co\int = 0; count < alien/2; -M-count) 

{ 

fprintf(fp, "tax %2i %4*lf\n", count / PZ, count % pz, 

*scale-M-) ; 

} 

fprintfCfp, "Channels C/D\n") ; 

for (count =0; count < aLen/2; -M-count) 

{ 

fprintf(fp, ••%2i %2i %4.1f\n", count / PZ, count % PZ, 

*scale-H-) ; 

} 

f close (fp) ; 

} 
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1, A method for illumination compensation for 
use in a system for reading surfaces having at least black 
and white areas, comprising the steps of: 

(a) storing reference reflected light 

5 intensity information corresponding to at least black and 
white areas in the form of a plurality of intensity val- 
ues, each one of said intensity values having an address 
defined by a value of at least one variable; 

(b) calculating from said reference re- 
10 fleeted light intensity information, a plurality of black 

threshold values and a plurality of white threshold val- 
ues, each one of said threshold values associated with 
selected values of each of said at least one variable ; 

(c) receiving operational reflected light 
15 intensity information in the form of a plurality of opera- 
tional reflected light intensity values, each one of said 
operational reflected light intensity values having an 
address defined by a value of each of said at least one 
variable ; 

(d) remapping each of said operational re- 
flected light intensity values to one of a plurality of 
output compensated reflected light intensity values, 
including the steps of remapping any of said operational ) 
reflected light intensity values equal to ott less than 
said black threshold value for its address to a black 
level output and remapping any of said operational re- 
flected light intensity values equal to or greater than 
said white threshold value for its address to a white 

level output; and 
30 (e) producing an output signal comprising 

said output compensated reflected light intensity values. 

2. A method as recited in claim 1, wherein 
said step of remapping further includes remapping any of 
said operational light intensity values less than said 
white. threshold value for its address and greater then 

. said black threshold value for its address to an output 
5 value equal to the absolute value of the difference be- 



20 



25 
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tween the white output level and the black output level, 
multiplied by the guotienr of said operational light 
intensity value less the black level threshold divided by 
the absolute value of the difference between the black 
10 level threshold and the white level threshold. 

3. A method as recited in claim 1, wherein 
said black output level is zero, and said white output 
level is a meiximum output level. 

4. A method as recited in claim 1, wherein 
said step of storing comprises storing reference reflected 
light intensity information corresponding to at least 
black and white areas in a plurality of profiles, each one 

5 of said profiles displaying said reference reflected light 
intensity information plotted against one of said vari- 
ables, each other of said variables being a selected 
constant for each profile. 

5. A method for illumination compensation in a 
system for processing of signals representing reflected 
light intensities from a surface, comprising the steps of: 

(a) receiving reflected light intensity 

5 information in the form of a plurality of intensity val- 
ues; 

(b) receiving, associated with each one of 
said reflected light intensity values, a value of at least 
one variable; 

10 (c) remapping each of said reflected light 

intensity values to one of a plurality of compensated 
reflected light intensity values, including the steps of 
remapping any of said reflected light intensity values 
equal to or less than a preselected threshold minimum 

15 reflected light intensity value associated with the value 
of each of said variables for said reflected light inten- 
sity value to a minimum level output and remapping any of 
said reflected light intensity values equal to or greater 
than a preselected threshold maximum reflected light 

20 intensity value, associated with the value of each of said 
variables for each reflected light intensity value, to a 
maximum level output. 

-41- 



wo 93/04442 PCT/L!S92/07054 

6. A method as recited in claim 5, wherein said 
srep of remapping further includes remapping any of said 
reflected light intensity values greater than said prese- 
lected threshold minimum reflective light intensity value 

5 associated with the value of each of said variables for 
each reflective light intensity value, and less than said 
preselected threshold maximum reflective light intensity 
value associated with the value of each of said variables 
for each reflected light intensity value, to an output 

10 value equal to the absolute value of the difference be- 
tween the maximum output level and the minimum output 
level, multiplied by the quotient of the reflective light 
intensity value less the threshold minimum light intensity 
value divided by the absolute value of the difference 

15 between the preselected threshold maximum reflective light 
intensity value and the preselected threshold minimum 
light intensity value. 

7. A method as recited in claim 6, wherein said 
minimum level output is zero and said maximum level output 
is a maximum possible output level. 

8. An apparatus for illumination compensation 
in a system for processing of signals representing re- 
flected light intensities from a surface, comprising: 

(a) means for receiving reflected light 

5 intensity information in the form of a plurality of inten- 
sity values; 

(b) means for receiving, associated with 
each one of said reflected light intensity values, a value 
of at least one variable; 

(c) memory means for remapping each of said 
reflected light intensity values to one of a plurality of 
compensated reflected light intensity values, including 
memory means for remapping any of said reflected light 
intensity values equal to or less than a preselected 
threshold minimum reflected light intensity value associ- 
ated with the value of each of said variables for said 
reflected light intensity value to a minimum level output 
and memory means for remapping any of said reflected light 
intensity values equal to or greater than a preselected 



10 



15 
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20 threshold maximxam reflected light intensity value associ- 
ated with the value of each of said variables for each 
reflected light intensity value, to a maximum level out- 
put. 

9. An apparatus as recited in claim 8, further 
comprising memory means for remapping any of said reflect- 
ed light intensity values less than said preselected 
threshold maximum light intensity value associated with 
5 the value of said variables for said reflected light 

intensity value and greater than said preselected thres- 
hold minimum. reflected light intensity value corresponding 
to a value of each of said variables for said reflected 
light intensity value, linearly to a value greater than 
10 said minimum level output and less than said maximum level 
output . 
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